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I.  Introduction 


The  project  entitled  "Study  of  Mutual  Coupling  between  Two  Slots  on  a 


Cylinder"  was  funded  by  Naval  Air  Systems  Command  under  Contract  N00019-76- 


M-0622  for  a period  of  four  months  (15  July  1976  to  15  November  1976).  The 


contract  monitor  was  Mr.  J.  Willis 


In  this  final  report,  we  summarize  (i)  personnel  (Section  II) 


(ii)  technical  results  (Section  III),  (iii)  plan  for  future  study  (Section  IV) 


and  (iv)  financial  statement  (Section  V) 


S.  W.  Lee,  Professor  of  Electrical  Engineering 


R.  Mittra,  Professor  of  Electrical  Engineering 


Rahmat-Samii 


Safavi-Naini,  Graduate  Research  Assistant 


(i)  UI  modal  solution:  We  have  derived  an  alternative  expression  for 


mutual  admittance  Y17  of  two  circumferential  or  axial  slots  on  a cylinder 


from  the  original  expression  given  by  Golden,  Stewart,  and  Pridmore-Brown 


IEEE  Trans. . vol.  AP-22,  pp.  43-48,  1974].  Advantages  of  this  UI  modal 


solution  are  as  follows:  (a)  Thre  is  no  need  to  introduce  a small  loss  in 


the  medium  for  the  purpose  of  circumventing  the  nonintegrable  singularity 


integrand  decays  exponentially  instead  of  being  highly  oscillatory  as  in 


the  GSB  expression.  Therefore,  as  long  as  z,  f 0,  the  UI  modal  solution  is 


more  suitable  for  numerical  calculations.  Detailed  discussion  of  the  UI 


modal  solution  is  given  in  Attachment  A 


2 


(ii)  UI  asymptotic  solution:  The  calculation  of  Y^  relies  on  a 

->  -> 

Green's  function  G(r,r')  for  the  surface  magnetic  field  due  to  a magnetic 
dipole  on  a cylinder.  We  have  developed  an  approximate  asymptotic  solution 
for  £(r,r'  ) valid  when  the  radius  of  the  cylinder  is  large.  A feature  of 
this  solution  is  that  G depends  not  only  on  the  longitudinal  curvature  along 
the  surface  ray,  but  also  on  the  surface  curvature  in  the  transverse  direction 


of  the  ray.  Using  our  solution  of  G,  for  circumferential  slots  has  been 
calculated  for  some  testing  cases.  The  results  are  in  excellent  agreement 
with  Y^'8  calculated  from  the  (exact)  modal  solutions  of  Hughes  or  UI. 

See  Attachment  B for  detailed  discussion. 


Plan  for  future  stud\ 


(i)  Using  the  UI  asymptotic  solution  for  £,  we  will  generate  more  data 
on  Y^2  80  as  to  ascertain  its  accuracy.  In  particular,  we  will  concentrate 
on  the  new  slot,  i.e..  Slot  B (a  = 0.5A,  b = 0.01A)  proposed  by  R.  C.  Hansen 
in  October  1976. 

(ii)  So  far  we  have  studied  only  the  coupling  between  circumferential 

slots,  which  depends  on  H, . Next  we  plan  to  investigate  the  coupling  between 

? 

the  axial  slots,  and  therefore  the  accuracy  of  the  H component  of  the  UI 

z 

asymptotic  solution.  It  should  be  noted  that  the  coupling  on  a cone 

depends  on  both  H.  and  H . 

<p  z 

(iii)  The  UI  asymptotic  solution  of  a cylinder  will  be  extended  to 
that  of  a cone.  This  can  be  done  only  after  we  gain  a better  understanding 
of  the  transverse-curvature  term  in  the  UI  solution.  We  do  not  expect  the 
latter  problem  to  be  an  easy  one. 


T.- 


Attachment  A 


CALCULATION  OF  MUTUAL  ADMITTANCE 
BETWEEN  TWO  SLOTS  ON  A CYLINDER* 


S.  Safavi-Naini  and  S.  W.  Lee 
Department  of  Electrical  Engineering 
University  of  Illinois  at  Urbana-Champaign 
Urbana,  Illinois  61801 


ABSTRACT 

The  mutual  admittance  between  two  identical  slots  on  the  surface  of  a 
cylinder,  which  can  be  expressed  in  terms  of  cylindrical  modal  functions, 
was  derived  by  Stewart,  Golden,  and  Pridmore-Brown.  Here  we  present  an 
alternative  expression  which  is  more  suitable  for  numerical  evaluation.  Our 
numerical  results  reveal  a somewhat  surprising  phenomenon:  the  coupling 

between  two  circumferential  slots  displaced  along  the  z-axls  is  signifi- 
cantly stronger  than  that  between  the  corresponding  slots  on  a planar  surface. 


* 


This  work  was  supported  by  The  Office  of  Naval  Research  under  Contract 
N00014-75-C-0293. 
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Consider  two  identical  circumferential  or  axial  slots  on  the  surface  of 

an  infinitely  long  cylinder  (Fig.  1).  Under  the  usual  assumption  of  a cosine 

aperture  field  distribution  ("one-mode"  approximation),  the  mutual  admittance 

has  been  determined  in  [1],  [2]  in  terms  of  cylindrical  modal  functions. 

Recently,  there  has  been  a renewed  interest  in  the  calculation  of 

because  (i)  it  is  a central  step  in  the  design  of  a conformal  array  [3],  [4], 

and  (ii)  it  provides  a convenient  check  for  various  asymptotic  (GTD) 

solutions  of  the  surface  ray  [5]  - [8].  The  original  expression  of  as 

given  in  [1],  [2]  is  not  suitable  for  numerical  calculations  when  the  slot 

separation  along  the  z-direction  is  large.  In  this  note,  we  will  present  an 

alternative  expression  for  Y12  which  overcomes  this  difficulty. 

Let  us  first  consider  the  circumferential  slots  shown  in  Fig.  1.  Denote 

the  center-to-center  separation  of  the  two  slots  by  ($q»2q)»  the  dimension  of 

each  slot  by  (a,b),  and  the  radius  of  the  cylinder  by  R.  For  exp(+jwt)  time 

* 

convention,  the  mutual  admittance  Y ^ is  given  in  Eq.  (8)  of  [2]  , which  reads 
in  the  present  notation. 
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The  multiplication  factor  2 in  the  definition  of  <j>jj  in  [2]  is  a misprint 
and  should  be  removed. 
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We  note  that  G contains  a finite  integral  and  can  be  evaluated  in  a straight- 
forward manner  by  standard  numerical  integration  techniques.  The  imaginary 
part  of  is  given  by 

“ cos  m<|>r 

V 

where  the  Integration  contour  is  shown  in  Fig.  2 and 
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The  computation  of  B as  given  in  (13)  can  be  quite  laborious  because  (i)  the 
integration  with  respect  to  k is  of  infinite  range,  and  the  factor  cos  k z_ 
is  highly  oscillatory  for  large  k^z^,  (ii)  W(m,kz)  has  nonintegrable  singulari- 
ties of  opposite  sign  on  both  sides  of  k * k_, (iii)  W(m,k  ) decays  slowly 

z u z 

with  respect  to  m and  kz> 

To  circumvent  the  above  difficulties  in  evaluating  B,  we  adopt  a 
method  introduced  by  Duncan  [9]  in  the  study  of  cylindrical  antenna  problems. 
Let  us  rewrite  (4)  as 


B - Im<  l 


cos  m$Q 


|m»0 

where 


-j  | F(m,kz)  sin  kzzQ  dkz  + J FGn.k^e^2*0  dk 


F(m,k  ) - [R(m,k  ) + jW(m,k  ) ]i(i(m,k  ) . 

2 2 2 Z 


'(5) 
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The  imaginary  part  of  the  first  term  inside  the  bracket  of  (5)  is 
,kz)  sin  kzzQ  dk. 


»{-j  f F(m,l 

l Jci 


- -|o"«<». 


k2)'l'(m.kz)  sin  kzzQ  dkz  . (7) 


In  order  to  compute  the  imaginary  part  of  the  second  term  of  (5),  the  integra- 
tion contour  is  deformed  into  C2  (Fig.  2)  according  to  the  theory  of 
complex  variables.  This  manipulation  leads  to 


f Jkz*0 

Im  F(m,k  )e  dk  • Im 

Jr  z * Jr 

L1  C2 


(8) 


Make  the  change  of  variable  k » jn  in  (8).  Substitution  of  the  resultant 


equation  and  (7)  into  (5)  gives 


Our  final  expression  for  Y.»  is  given  in  (2),  with  its  real  part  G in  (3)  and 


its  imaginary  part  B in  (9) 


but  also  B is  determined  by  R(m,k  ),  which  is  much  simpler  than  W(m,k  ) defined 


in  (4b).  (ii)  G contains  only  a finite  integral.  (iii)  The  infinite  inte 


gral  in  B,  i.e.,  the  second  integral  in  (9a),  contains  an  exponentially 


decaying  factor  exp[-(z 


faster  convergence  is  the  evaluation  of  B.  This  is  in  contrast  to  the 


original  expression  of  given  in  (1).  (iv)  There  is  no  nonintegrable 


singularity  in  (3)  or  (9) 


The  same  method  applies  to  the  derivation  of  an  alternative  expression 
_ for  two  axial  slots  (a  < b as  shown  in  Fig.  lb) . We  give  below  only 


the  final  result 


♦ ■ 2 sin  (a/2R) 
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Some  typical  numerical  data  of  Y^  calculated  from  our  expressions  in 
(2)  and  (10)  are  presented  in  Tables  I and  II.  The  parameters  are  f » 9 GHz, 
R * 1.991",  a * 0.9"  and  b *»  0.4"  for  circumferential  slots,  and  a = 0.4" 
and  b ■ 0.9"  for  axial  slots.  Y ^ is  listed  in  (db,  phase  in  degree),  where 
db  ■ 20  l°g^Q  with  Y^  in  mho.  When  Zq  is  small  (less  than  4"),  our 

results  are  in  excellent  agreement  with  those  given  in  [1],  [2],  [4].  For 
large  Zq,  the  original  expressions  of  Y^  [1]»  [2]  become  very  slowly  con- 
vergent and  no  numerical  results  have  been  reported. 

With  our  new  expressions,  we  have  done  extensive  numerical  studies  of 

the  asymptotic  behavior  of  Y^  for  large  Some  results  are  presented 

in  Figs.  3 and  4,  where  Iy^I  of  circumferential  and  axial  slots  is  plotted 

as  a function  of  z^  for  ^ * 0 and  a cylinder  of  radius  R = 1.991".  We 

observe  that  (i)  Iy^I  of  axial  slots  (Fig.  3)  on  a cylinder  is  very  close 

to  that  on  a plane  (R  -►  ®) , and  (ii)  lY^I  of  circumferential  slots  (Fig.  4) 

-1/2  -1 

on  a cylinder  asymptotically  decays  as  (^qZq)  in  contrast  to  the  (k^z^) 
behavior  on  a plane.  The  significantly  stronger  coupling  between  circum- 
ferential slots  on  a cylindrical  surface  as  compared  to  that  on  a planar 
one  is  rather  unexpected,  and  it  has  a very  important  implication  on  the 
current  GTD  formulas  for  calculating  surface  rays.  The  latter  aspect  is 
studied  in  [8]. 
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ABSTRACT 

A simple,  approximate  expression  for  the  surface  magnetic  field  due  to 

a magnetic  dipole  on  a conducting  circular  cylinder  is  obtained.  This 

solution  is  asymptotic  for  a large  cylinder  radius,  and  is  uniformly  valid 

everywhere  on  the  cylindrical  surface  including  the  penumbra  and  the  deep 

shadow.  In  the  limit  that  the  cylinder  radius  is  infinite,  it  becomes 

identical  to  the  known  exact  solution  of  a dipole  on  a conducting  plane. 

For  a surface  ray  propagating  in  parallel  to  the  axis  of  the  cylinder,  the 

-1/2 

transverse  surface  magnetic  field  is  found  to  vary  asymptotically  as  (ks)  , 
where  s is  the  distance  from  the  source.  This  behavior  is  distinctively 
different  from  the  (ks)  * variation  of  the  surface  ray  on  a plane,  and  is 
explained  in  terms  of  the  dependence  of  the  surface  curvature  in  the  binormal 
direction  of  the  ray.  We  apply  our  solution  to  the  mutual  coupling  problem 
between  two  slots  on  a cylinder,  and  obtain  results  which  are  in  excellent 
agreement  with  those  calculated  from  the  exact  modal  solution.  A comparison 
of  the  present  solution  with  two  other  asymptotic  (GTD)  solutions  is  also  given. 
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1.  INTRODUCTION 


This  paper  considers  a high-frequency  diffraction  problem  by  a perfectly 
conducting  cylinder  as  sketched  in  Fig.  la.  For  a given  magnetic  dipole  at 
Q'  on  the  surface  of  the  cylinder,  the  problem  is  to  find  the  surface  magnetic 
field  (surface  current  density)  everywhere  when  kR  is  large  (R  is  the  radius 
of  the  cylinder  and  k » 2n/A).  The  motivations  for  our  study  are  the  following 
two:  (i)  The  solution  of  the  cylinder  problem  constitutes  a central  step  in 

calculating  the  mutual  coupling  between  two  slots  on  the  surface  of  a 
cylinder  [l]-[8].  (ii)  More  importantly,  the  cylinder  problem  is  a so-called 
"canonical  problem"  in  GTD  [9]— [11] . Once  its  solution  is  known,  it  may  be 
generalized,  by  following  the  recipe  of  the  GTD,  to  give  the  asymptotic  solu- 
tion of  the  surface  magnetic  field  on  any  convex,  smooth,  conducting  surface. 

The  cylinder  problem  has  an  exact  modal  solution,  which  is  in  the  form 
of  an  infinite  series  with  each  term  containing  an  infinite  integral  [4],  [5]. 
For  a large  kR,  this  solution  is  very  slowly  convergent  and  becomes  less 
useful.  Two  asymptotic  solutions  exist  in  the  literature:  one  given  by 

Hwang,  Kouyoumjian,  and  Pathak  [5],  (6]  (hereafter  referred  to  as  the  OSU 
solution),  and  the  other  by  Chang,  Felsen,  Hessel  and  Shmoys  (3],  [4]  (the 
PINY  solution).  Both  are  approximately  deduced  from  the  exact  modal  solution 
under  the  condition  kR  * ■».  In  the  present  paper,  we  offer  a third  asymptotic 
solution  (the  UI  solution),  which  gives  the  surface  magnetic  field  everywhere 
from  the  source  point  to  the  deep  shadow  in  a single  expression,  and  is  based 
on  a classical  work  by  Fock  in  1949  (Chapter  12  of  [12]). 

The  organization  of  this  paper  is  as  follows.  In  Sections  2 through  4, 
the  final  form  of  the  UI  solution  is  stated  and  compared  with  those  of  OSU 
and  PINY.  In  Section  3,  the  three  asymptotic  solutions  are  applied  to  the 


evaluation  of  the  mutual  admittance  between  two  slots  on  a cylinder.  Their 


results  are  compared  with  those  calculated  from  the  exact  modal  solutions  [2], 
17),  [ 8 ] . Section  6 describes  the  derivation  of  the  I'l  solution.  Finally, 
a conclusion  is  given  in  Section  7.  Some  formulas  of  Fock  functions  used  in 
the  text  are  listed  in  the  Appendix. 
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2.  UI  SOLUTION  FOR  THE  SURFACE  MAGNETIC  FIELD 

At  point Q'  on  the  surface  of  the  cylinder  (Fig.  la),  there  is  a tangential 
magnetic  dipole  source  described  by  a magnetic  current  density  (for  exp  +jwt 
time  convention) 

hr)  « M | 6(r  - R)6(4>)6(z)  (2.1) 

where  M is  the  magnetic  dipole  moment  (M  • r = 0),  and  (r  = R,4>  = 0,z  = 0) 
are  the  cylindrical  coordinates  of  Q'.  The  problem  is  to  determine  H at 
another  point  Q = (R,ij>,z)  on  the  surface  under  the  assumption  that  kR  is  large. 

First,  let  us  introduce  several  parameters.  According  to  GTD  [9],  [11], 
the  dominant  contribution  of  H at  Q is  the  field  on  the  surface  ray  from  Q'  to  Q. 
The  surface  ray  is  a geodesic  on  the  conducting  surface,  and  in  the  present 
case  is  a helical  path  (Figure  1).  The  arclength  of  the  surface  ray  is 


s = 


(2.2) 


The  tangent,  normal,  and  binormal  of  the  surface  ray  are  (t'.-n'.-b')  at  Q', 

A A A A A A 

and  (t,-n,-b)  at  Q.  Thus,  (t,n,b)  form  a moving  trihedron  iilong  a surface  ray, 
pointing  toward  the  longitudinal  and  two  transverse  directions.  At  any  point 
on  the  surface  ray,  the  curvature  of  the  conducting  surface  is  described  by 
two  parameters: 

Rj.  » the  radius  of  curvature  in  the  direction  of  t (or  that  in  the 
longitudinal  direction  of  the  surface  ray),  and 


R^  = the  radius  of  curvature  in  the  direction  of  b (or  that  in  the 


transverse  direction  of  the  surface  ray). 
On  a convex  surface,  both  Rt  and  R^  are  nonnegative, 
a conducting  cylinder,  one  has 
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2 . 
cos  6 


*b  ' 


sin2  6 


(2.3) 


where  0 is  measured  from  R4>-axis  in  Fig.  lb,  and  takes  a value  between  0 and 
2 ir.  The  large  parameter  for  our  asymptotic  expansion  is 

1/3 


■j  kR  ) 

x2  t) 


(2.4) 


Thus,  our  solution  is  an  approximate  (not  rigorous)  asymptotic  solution  valid 
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for  m > up  to  and  including  terms  of  0(m  ).  Furthermore,  let  us  introduce 

a distance  parameter 


/ 0\ 1/3 

£ - jp  - k/2R2*  s 

Kt  v c / 


ks 

2m2 


(2.5) 


which  is  the  arclength  normalized  by  k and  R^ . Note  that  £ * 0 defines  the 
lit  region  (0  = tt/2),  £ < 1 defines  the  penumbra  region,  and  C >>  1 defines 
the  deep  shadow.  Our  solution  is  uniformly  valid  for  all  £ » 0. 

Due  to  the  point  source  in  (2.1),  our  final  asymptotic  solution  for  the 
magnetic  field  on  the  surface  is  given  by 


-4  „ 


H(Q)  = M • (b'bl^  + t'tHt)  + 0(m  ),  m * * 

where  the  transverse  component  is 

- 2 

Hb(Q)  % f1  “ ii)v(^  " JE)  u(°  + j(^  kRt)~2/3v'(£) 

+ j(/2  kRt)'2/3(Rt/Rb)u'(0  G(s)  , 

- 

the  longitudinal  component  is 


(2.6a) 


(2.6b) 


Ht(Q)  % (i)  v<°  + I1  ~ kilu(0  + j(/2  kv~2/3l,,(° 


G(s)  , (2.6c) 


and  the  function  G(s)  is 


(120n) 


v and  u are  defined  in  the  Appendix,  and  v 


is  the  derivative  of  v.  We  emphasize  that  (2.6)  is  an  approximate  solution 


For  one  thing,  we  are  not  even  able  to  show  that  (2.6a)  contains  all  the  terms 


up  to  0(m  ),  let  alone  other  more  subtle  questions 


Let  us  consider  several  limiting  cases  of  the  UI  solution  given  in 


(2.6).  If  the  radius  of  the  cylinder  becomes  infinite 


the  use  of  (A-12)  through  (A-16)  in  the  Appendix  in  (2.6)  leads  to 


When  (2.8)  is  substituted  into  (2.6a),  we  find  that  H in  (2.6a)  is  identical 


to  the  exact  solution  of  the  surface  field  due  to  a magnetic  dipole  on  a flat 


ground  plane  [A],  [13] 


The  second  limiting  case  occurs  when 


We  find  from  (2.6)  that  H is  again  given  by  (2.8b)  but  H becomes 


In  terms  of  the  planar  solution  in  (2.8a),  we  may  rewrite  (2.10a)  as 


planar  8 /2ir 


The  result  in  (2.10)  is  most  interesting,  and  in  fact  somewhat  surprising. 
The  surface  ray  traveling  in  direction  0 = rr/2  (Fig.  1)  is  a straight  line 
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(kRt  ■+  '») . However,  due  to  the  finite  curvature  in  the  binormal  direction 
(R^  ■ R),  on  the  cylindrical  surface  differs  from  its  counterpart  on  a 
planar  surface  by  the  additional  term  in  (2.10b).  At  a large  distance  away 
from  the  source  (ks  -*•  <*>)  in  the  direction  0 = n/2,  and  for  a fixed  kR,  we 
find  that  on  a planar  surface  and  that  on  a cylindrical  surface  are  given 
by,  respectively. 


|Hb(Q)] 


-jks 


planar 


“b«> 


A 


^ B 


ks 


1 e 


-jks 


-jks 


kR 


/ks 


+ A 


ks 


(2.11) 

(2.12) 


where  A and  B are  constants  independent  of  s and  R.  Thus,  for  large  ks,  H. 
on  a cylinder  is  stronger  than  that  on  a plane.  Such  a phenomenon  was  first 
reported  by  Hasserjian  and  Ishimaru  [14],  and  later  by  the  authors  of  [4] 
and  [7].  Those  previous  workers,  however,  have  not  explained  the  phenomenon 
in  terms  of  R^,  the  radius  of  curvature  in  the  binomial  direction,  as  we  did 
in  (2.6b). 

As  a third  limiting  case,  let 

£ > ® (2.13) 

which  occurs  when  observation  point  Q is  in  the  deep  shadow.  Making  use  of 
(A-7)  through  (A-ll),  we  have  from  (2.6)  that 


.2  2/3  . 

„ /„*  . k cos  6 

Hb(Q)  1/3 — Til  exp 

1528(kR)i/J(ks)i/z 


0 . 88£,  - j (yj  + 0.5U  + ksj 


f -►  IX> 


Vq)  ' ft  VQ)  • 5 


(2.14a) 


(2.14b) 


Therefore,  in  the  deep  shadow,  the  field  is  a slow  wave  and  decays  exponentially 
along  the  surface  ray. 


/ 
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In  applying  the  formulas  in  (2.6)  to  the  mutual  admittance  calculation* 
[l]-[8],  explicit  field  expressions  are  needed  for  axial  and  circumferential 
dipoles.  They  are  listed  below. 

Axial  dipole:  M - z 


Hz(Q)  - ^ cos2  6 + Ht  sin2  0 


~ (v(4) 


{• 

+ j(/2  kRt) 


2 i 

cos  0 - cos  20 


r 


-2/3 


v'  (4)  cos  0 + 1 


+ i)  U’(C)  811,2  6]} 


Circumferential  dipole:  M 


(2.15a) 

+ (Al 

2 

cos  0 

G(s) 

(2.15b) 

(2.16a) 

H^(Q)  * sin2 0 + Ht  cos2  0 


~ ^v(C) 

+ j(/2  kR  )_2/3 


sin  6 + ^ cos 


29j * (i) U<u  cos2  9 i1  • li) + lii) s1"2  9 


v'(4)  sin2  0 + ^tan^  0 


2 

cos  0 


) 


G(s)  . 
(2.16b) 


In  the  limiting  case  4+0  (either  kR  + 00  or  0 + n/2),  (2.15b)  becomes 


Hz(Q)  * 


cos2  0 (2  - 3 cos2  0)  (l  - j^) 


G(s) , 4+0  (2.17) 


and  (2.16b)  becomes 


H^(Q)  % [sin2  0 + £ (2-3  sin2  0)  I1  ■ i)  + wJg(s),  4 + 0 . (2.18a) 


where 


, if  kR  + 00 


(2.18b) 


W 


, lf  e .,/2 


V 


(2.18c) 


Note  that  in  the  limit  kR  + <*>,  (2.17)  and  (2.18)  recover  the  exact  solutions 
for  dipoles  on  a planar  surface. 


1 
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3.  COMPARISON  WITH  OSU  SOLUTION 

For  the  same  problem  of  a magnetic  dipole  on  a cylinder,  an  asymptotic 

(GTD)  solution  was  derived  from  canonical  problems  at  OSU.  That  solution, 

given  in  Section  D of  [5]  or  Eqs.  (6.64)  and  (6.80)  of  [11],  can  be  also 

* 

written  in  the  form  of  (2.6a)  with 


Hb(Q)  * v(£)G(s) 

Ht(Q)  % (i|u(OC(s)  . 


(3.1a) 

(3.1b) 


For  the  extreme  case  £ -*■  0 (either  kR  -*•  °°  or  0 -*■  tt / 2 ) , the  use  of  (2.7)  in 
(3.1)  leads  to 


^(Q)  G(s) , £ 0 

Ht(Q)  * (i)G(s),  5.0  . 


(3.2a) 

(3.2b) 


For  the  other  extreme  case,  £-*•««,  the  use  of  (A-7)  and  (A-8)  in  (3.1)  leads 
to  an  identical  to  that  in  (2.11a),  and 

2 1/2  2 

v»  ' li)  -- - - j If  ■ * i-iw  * HD-  * * - • 

(3.3) 

The  OSU  solutions  in  (3.1),  (3.2),  and  (3.3)  should  be  compared  with  our  so- 
lutions in  (2.6),  (2.8),  (2.10)  and  (2.14).  Several  remarks  are  in  order. 

(i)  In  the  limit  kR  •+  ®,  our  solution  in  (2.8)  is  identical  to  the  known 

exact  solution.  On  the  other  hand,  in  (3.2a)  recovers  only  the  term  of 

-1  -2  -3 

(ks)  , but  not  the  terms  of  (ks)  or  (ks)  . The  latter  terms  are  important 

for  the  field  near  the  source.  For  Ht  in  (3.2b),  a factor  2 is  missing  in 
_2 

the  term  of  (ks)  , as  pointed  out  in  [5]. 


Note  the  corresponding  notations  in  [5]  and  here:  ip(£)  -*■  2e  ^^v(£). 


¥<0 


e 3/2^^^  a -►  0,  a -*•  R,  p R , and  t -*■  s. 

8 • 


f 
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(11)  For  a large  but  finite  kR,  the  two  solutions  in  (2.6)  and  (3.1)  do 
not  agree.  In  the  penumbra  region  (£  < 1),  the  UI  solution  in  (2.6)  should 
be  more  accurate  because  of  (i).  In  the  deep  shadow,  both  H^s  are  given  by 
(2.14a),  but  the  two  solutions  for  Ht  in  (2.14b)  and  (3.3)  are  completely 
different.  In  Sections  4 and  6 we  will  show  that  (2.14b),  not  (3.3), 
agrees  with  the  PINY  solution,  and  gives  more  accurate  numerical  results. 

(iii)  Because  of  the  observations  in  (i)  and  (ii),  it  appears  that  Hf 
given  in  (3.1b)  is  not  accurate. 

(iv)  For  a fixed  kR  and  in  the  direction  9 = tt/2,  in  (3.1a)  becomes 
asymptotic  for  a large  ks. 


V<» 


which  should  be  compared  with  the  UI  solution  in  (2.12).  We  note  that  the 
term,  attributed  to  the  curvature  in  the  binormal  direction  of  the  surface 
ray,  is  absent  in  (3.4). 

(v)  For  acoustic  diffraction  by  a cylinder,  the  functions  (v,u)  arise 
when  the  boundary  condition  is  (hard,  soft).  We  note  from  the  OSU  solution 
in  (3.1)  that  H^  depends  on  the  "hard"  function  v,  while  depends  on  the 
"soft"  function  u.  Such  a separation,  however,  is  not  possible  for  the  UI 
solution  in  (2.6). 

(vi)  In  Section  C of  [5],  Hwang  and  Kouyoumjian  modified  their  solution 


of  Ht  in  (3.1b)  to  read 


Ht(Q)  * T(i)u^>G(s)  ' 


Here  the  additional  factor  T for  the  present  cylinder  problem  is 


where  9'  = sin 


-1  M • z 


i 


(3.6) 


They  attributed  the  presence  of  T to  the  torsion  of  the  surface  ray.  However, 
there  exist  several  apparent  difficulties  in  connection  with  T:  (a)  For  a 

nonaxial  dipole  (9'  f n/2),  T becomes  infinite  as  0 + tt/2.  (b)  If  T is  indeed 

a torsion  factor,  it  should  be  reduced  to  unity  for  a torsionless  ray  propa- 
gating in  the  $ direction  on  a cylinder  (0  = 0).  However,  T calculated  from 
(3.6)  does  not.  (c)  An  arbitrarily  oriented  dipole  on  the  surface  of  a 
cylinder  may  be  resolved  into  an  axial  dipole  and  a circumferential  dipole. 

We  may  calculate  fields  due  to  each  dipole  separately,  and  later  superimpose 
them  for  the  original  solution.  However,  if  formulas  in  (3.5)  and  (3.1a)  are 
used,  such  a superposition  procedure  does  not  recover  the  original  solution. 
Because  of  the  above  difficulties,  we  will  use  (3.1b),  not  (3.5),  for  all 
the  subsequent  discussion  of  the  OSU  solution.  (In  [3],  [4]  it  is  (3.5), 
not  (3.1b),  that  was  used  in  all  the  numerical  calculations.) 

We  would  like  to  emphasize  that  the  OSU  solution  represents  one  of  the 
very  first  efforts  to  apply  the  ray  technique  to  the  surface-field  calculation. 
Their  solution,  while  inadequate  in  some  situations,  has  produced  many  useful 
results  [10],  [11],  and  more  importantly,  has  laid  a conceptual  framework 
from  which  the  more  refined  works,  e.g.,  PINY  and  UI  solutions,  are  deduced. 
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4.  COMPARISON  WITH  PINY  SOLUTION 


Another  asymptotic  solution  (PINY  solution)  for  the  cylinder  problem  was 
deduced  from  the  exact  modal  solution  [3],  [4],  and  is  given  by 


Axial  dipole:  M = z 


Hz(Q)  -v  v(4)  cos2  8 + (2  - 3 cos2  6)  G(s) 

- 


Circumferential  dipole:  M = <J> 


MQ>  * 


f r — 

sv(4)  sin2  0 + (1-3  sin  0) 

' L 

+ sec2  0 [u(U  - sin2  0 V1(C) A G(s) 


(4.1) 


(4.2) 


where  v^(4)  is  defined  in  the  Appendix.  In  the  limiting  case  4 -►  0 (either 
kR  -*■  °°  or  0 -+  n/2),  (4.1)  becomes 


Hz(Q)  -v  cos2  0 + ^ (2  - 3 cos2  0)  G(s),  4 -►  0 (4.3) 


and  (4.2)  becomes 


where 


H (Q)  -v  sin2  0 + -jj:  (2  - 3 sin2  0)  + W G(s),  4 + 0 


, if  kR  -*■  00 


1/2  ,,  ,,  % 1/2 
(f,  e-J3lf/4  , if  0 + w/2 


(4.4a) 


(4.4b) 


(4.4c) 


The  PINY  solutions  in  (4.1)  through  (4.4)  should  be  compared  with  the  UI 
solution  in  (2.15)  through  (2.18).  The  following  observations  are  made: 


There  are  several  slightly  different  formulas  given  in  (4].  The  ones  pre- 
sented here  are  the  "full  formulas"  taken  from  Eqs.  (101)  and  (111)  of  [4]. 
Note  the  corresponding  notations  in  [4]  and  here:  D -»  s,  x -*■  4,  -*•  v, 

v^  -*■  v^,  and  Uq  -*■  u. 
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(i)  In  the  limit  kR  -*■  °°,  the  PINY  solution  recovers  the  exact  solution 

-1  -2  -3 

for  a planar  surface  in  terms  of  (ks)  and  (ks)  , but  not  in  terms  of  (ks) 

(ii)  For  a fixed  kR  and  in  the  direction  0 = tt/2,  both  H of  the  PINY 

solution  in  (4.4)  and  that  of  the  UI  solution  in  (2.18)  vary  asymptotically 

-1/2  -1 
as  (ks)  , which  is  distinctively  different  from  the  (ks)  behavior  in  the 

planar  solution.  However,  the  factor  W in  (2.18c)  contains  an  extra  factor 

(3/4)  when  compared  with  its  counterpart  in  (4.4c). 

(iii)  In  the  deep  shadow  C 00 , both  H and  H,  of  the  PINY  solution  decay 

z <p 

exponentially  according  to  the  "hard"  function  v(£),  which  is  in  agreement 
with  our  solution  in  (2.11).  In  particular,  at  0 = 0,  H.  in  (4.2)  of  the 
PINY  solution  becomes 

VQ)  ’ HtW) '•  |i)  "!.«»-  6 ■ 0 (4-5> 

where  H^(Q)  is  given  in  (2.14a).  Note  that  (4.5)  agrees  with  the  UI  solution 
in  (2.14b),  but  disagrees  with  the  OSU  solution  in  (3.3).  It  has  been  demon- 
strated in  [4]  (see  Figures  III-2  and  III-13)  that  the  use  of  (4.5)  for  mu- 
tual admittance  calculations  between  slots  shows  good  agreement  with  an 
exact  numerical  solution,  while  the  result  calculated  from  (3.3)  deviates 
markedly  from  the  exact  solution. 

(iv)  Hz  in  (4.1)  depends  on  the  "hard"  function  v(£)  only,  whereas  our 
solution  in  (2.15)  depends  on  both  v(C)  and  u(£) . 

(v)  In  Eq.  (123)  of  [4],  an  expression  for  H was  given  for  an  arbitrarily 
oriented  dipole.  It  has  the  form 

H(Q)  ^ M • [b'bA  + t'tB  + $'$C]  (4.6) 

which  should  be  compared  with  (2.6a).  The  presence  of  the  cross  term  C in 
(4.6)  is  unique. 
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Next,  we  present  a numerical  comparison  of  H,  for  all  three  solutions: 

9 

the  UI  solution  in  (2.16b),  the  OSU  solution  in  (3.1),  and  the  PINY 

solution  in  (4.2).  For  a cylinder  with  radius  kR  = 9.5325,  H,  on  three 

9 

surface  rays  (6  » 0°,  45°,  and  90°)  is  displayed  in  Figs.  2 to  5 as  a 

function  of  ks,  the  distance  from  the  source.  The  magnitude  of  H is 

9 

normalized  by  H,  (ks  = 0.2,6)  of  the  UI  solution,  whereas  the  normalized 
9 

i ks 

phase  is  equal  to  Arg(je  H^).  Discussions  of  these  numerical  results  are 
given  below:  (a)  As  ks  -*■  0,  only  the  UI  solution  exhibits  the  correct 

_3 

behavior  (ks)  . This  fact  explains  the  marked  disagreement  among  the  three 

solutions  in  the  range  0 < ks  < 1.  (b)  In  the  direction  6=0  (Fig.  2),  UI 

and  PINY  solutions  converge  to  each  other  for  ks  > 6 (or  £ > 1) , while  the 

OSU  solution  is  much  too  small,  as  noted  in  (iii)  above.  (c)  In  the  direction 

-1/2 

6 = 90°  in  Fig.  5,  both  the  UI  and  PINY  solutions  vary  as  (ks)  , while 
the  OSU  solution  varies  as  (ks)  ^ . (d)  Overall,  the  three  solutions  are 

significantly  different.  In  applications  where  accuracy  within  one  or  two 
db  is  required,  the  selection  of  a proper  formula  becomes  crucial. 

To  present  an  overall  view  of  the  surface  magnetic  field  due  to  a circum- 
ferential dipole  (M  = 9),  three-dimensional  plots  of  |h^|  and  | | calculated 
from  the  UI  solution  (2.6)  are  given  in  Figure  6.  Fields  in  both  plots  are 

normalized  by  I H , I at  ks  = 0.43  and  6 = u/2,  which  has  a value  0.00503  a/m 
9 

and  is  at  zero  db.  Field  values  above  7.475  db  are  not  shown  in  these  plots. 


15 


NORMALIZED  PHASE 


' OSU (PHASE) 


JJI  (PHASE) 


.PINY 


\OSU  MAG 

\ 

\ 

^ * s\ 


-17°  g 
z 


6 * 90° 
kR  = 9.5325 
OdB  = 8.386x10' 


• ^PXNY  (PHASE) 

s \ 
s \ 

V \ 


Figure  5.  along  a surface  ray  in  0 * 90°  direction  on  a cylinder  vs.  the 
distance  from  the  ^-directed  magnetic  dipole. 
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5.  APPLICATION:  MUTUAL  ADMITTANCE  OF  SLOTS 

Let  us  apply  the  asymptotic  solutions  given  in  the  previous  three  sections 
to  the  calculation  of  the  mutual  admittance  Y^  between  two  slots  on  a 
cylinder.  This  is  done  not  only  because  the  Y^  calculation  is  an  important 
practical  problem,  but  also  because  of  the  existence  of  the  exact  modal  solu- 
tion for  which  provides  a convenient  numerical  check  of  the  accuracy 

of  the  asymptotic  solutions. 

Referring  to  Fig.  7,  let  us  consider  two  identical  circumferential  slots 
on  the  surface  of  a cylinder.  Under  the  condition  that 

(a/A)  = 0.5  and  (b/A)  « 1 (5.1) 

it  is  reasonable  to  approximate  the  aperture  field  of  the  slot  1 (or  slot  2) 
by 

E1(y,z)  = cos  ~ y , for  |y|  < a/2,  \z\  < b/2  (5.2) 

where  y * R$  is  the  angular  distance  along  the  ^-direction,  and  is  the 
voltage  difference  across  the  slot.  Once  the  "one-mode"  approximation  in 
(5.2)  is  accepted,  it  can  be  shown  [2],  [4]  that  the  mutual  admittance  between 
the  two  slots  in  Fig.  7 is  given  by 
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^2 

ab 


dyl  dzl  dy2  dz2  COS  a yl  COS  a (y2  “ yo*  * 

J A * A _ ^ ~ v 


(5.3) 


Here  A^  and  are  the  apertures  of  the  slots,  and  (Yq»zq)  are  their  center- 
to-center  distances.  The  Green's  function  g in  (5.3)  represents  the  at 
(y 2*z2^  ^ue  to  3 unit-strength,  ^-directed  magnetic  dipole  at  (y^,z^).  For 
UI,  OSU,  and  PINY  asymptotic  solutions  of  Y^,  g is  equal  to  in  (2.16b), 
(3.1),  *nd  (4.2),  respectively. 
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In  addition  to  the  three  asymptotic  solutions  of  Y^.  there  exists  an 

expression  of  Y^  in  terms  of  cylindrical  eigenfunctions  as  given  in  Eq.  (8) 

* 

of  [2].  Within  the  "one-mode"  approximation,  this  expression  is  exact,  and 
contains  an  infinite  series  with  each  term  in  the  series  being  an  infinite 
integral.  This  expression  has  been  evaluated  by  two  different  numerical 
techniques.  We  call  the  one  reported  in  [2],  [7]  the  Hughes  modal  solution, 
and  the  one  in  [ 8 ] the  UI  modal  solutions. 

Now  we  will  present  some  numerical  results  calculated  from  the  two  (exact) 
modal  solutions,  and  the  three  (approximate)  asymptotic  solutions.  The 
parameters  of  the  cylinder  and  the  slots  are 

f = 9 GHz,  a = 0.9",  b = 0.4",  R = 1.991"  . (5.4) 

As  a function  of  slot  separation  (z  ,$q),  is  listed  in  Tables  I to  III  in 
(db  value,  phase  in  degree),  where  the  db  value  is  calculated  from  the 
relation 

db  = 20  log1()  |Y12/Y11|  , |YU|  = 1.7075  * 10'3  mho  . (5.5) 

From  the  comparison  made  in  Tables  I to  III,  we  conclude:  (a)  Considering 

the  numerical  integration  error  involved,  the  UI  asymptotic  solution  is  in 
excellent  agreement  with  the  (exact)  modal  solution.  (b)  The  PINY 
asymptotic  solution  is  reasonably  accurate,  while  OSU  is  not.  In  addition 
to  the  data  presented  in  this  paper,  we  have  calculated  nearly  100  different 
Y^'s  over  a wide  range  of  parameters.  The  above  two  conclusions  hold  for 
all  the  calculations  provided  that  kR  _>  5. 

As  discussed  in  Section  2,  the  surface  ray  propagating  in  the  direction 
0 = it/ 2 on  a cylinder  has  a stronger  H,  than  its  counterpart  on  a plane. 


There  is  a misprint  in  the  definition  of  <}>b . The  correct  one  should  read 
<|>k  = sin-^  (b/2pg).  It  should  be  also  pointed  out  that  the  corresponding 
definition  used  in  (5.2)  is  = (b/2pg)>  For  numerical  data  presented  in 
this  paper,  the  two  slightly  different  definitions  of  d>b  have  negligible 
effects. 
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Since  along  the  E-plane  ($q  * 0)  is  proportional  to  Hb,  this  phenomenon 
is  also  seen  in  Figure  8,  where  we  plot  the  ratio 

Y^  on  a cylinder  with  radius  R 
Yj^  00  3 plane 

as  a function  of  R for  = 8"  and  4>q  = 0.  We  note  that  the  convergence 
rate  of  the  cylindrical  Y^  to  the  planar  Y^  is  n°t  as  rapid  as  one  would 
normally  expect.  For  example,  at  kR  = 50,  the  cylindrical  Y^  Is  still 
about  10  percent  higher  than  the  planar  one. 


y on  a cylinder  as  a ^unc^^-on  the  rad^us  ^ the  cyiinder*  ^ 
normalized  by  Y,,  on  a plane  which  is  5.37  x 10  exp(j53.55°)  mho. 


I 
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6.  DERIVATION  OF  UI  SOLUTION 

The  UI  solution  given  in  (2.6)  will  now  be  derived.  The  derivation  is 
based  on  the  following  observation.  Consider  a surface  ray  (a  geodesic) 
on  an  arbitrary  convex  surface.  In  the  neighborhood  of  the  ray,  the  geo- 
metrical properties  of  the  surface  depend  on  two  radii  of  curvature:  Rfc  in 

the  direction  of  tangent  t and  in  the  direction  of  binormal  b.  It  is 
well-known  that  the  dominant  asymptotic  solution  of  a surface  field  depends 
on  Rt>  not  R^  (see  the  discussion  in  pp.  192-193  of  [11]).  Thus,  for  the 
purpose  of  determining  the  dominant  term  of  the  surface  field  in  the  direction 
of  0 in  Fig.  1,  the  cylinder  may  be  replaced  by,  for  example,  a sphere  of 
radius  Rt*  In  other  words,  our  problem  is  then  to  solve  the  radiation  of  a 
magnetic  dipole  on  a perfectly  conducting  sphere. 

In  studying  wave  propagation  around  the  earth  in  1949,  Fock  considered 
the  radiation  of  a horizontal  electric  dipole  on  a lossy  dielectric  sphere  [12]. 
He  first  obtained  an  exact  solution  in  terms  of  spherical  harmonics,  and 
next  extracted  the  dominant  high-frequency  terms  from  it.  By  the  duality 
principle,  Fock's  solution  can  be  converted  to  the  one  that  is  sought  by 
us.  We  will  give  below  the  solution  after  the  conversion  is  done. 

Consider  a magnetic  dipole  in  the  x-direction  (M  = x),  located  on  a 
perfectly  conducting  sphere  of  radius  Rt  (Fig.  9).  The  total  field  can  be 
derived  from  a magnetic  Hertz  potential  U and  an  electric  Hertz  potential  V. 

The  asymptotic  solutions  of  U and  V are  found  to  be  [from  (7.09)  and  (7.10) 
on  p.  252  of  [12]]* 


The  divergence  factor  DF  of  the  surface  ray  on  a sphere  has  been  removed  in 
(6.1)  [DF  » (0/sin  0)^^].  Note  the  corresponding  notations  in  [12]  and 
here:  i -+  (-j),  u -*•  U,  v -*•  V,  a -*•  Rf,  and  x ■*  £. 
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U<?)  * G(S)  -^-*-a-J*/4(£/.)1/2 

2k^Rtm 

• |r  e_jCt(2)w2(t)[wi(t  - *i> 


(t) 


(t)  2 ' ’1 


w (t  - y, ) dt  (6.1a) 


Hn 


* G(s)  ^-A  V(C) 
Rt  k Rt 


(6.1b) 


where  = (k/m)(r  - R^),  and  (r,0,4>)  are  the  spherical  coordinates  of  r. 
Parameters  G,  m,  and  £ in  (6.1)  are  defined  in  Section  2,  while  the  others 
are  in  the  Appendix.  The  magnetic  field  is  calculated  from  the  relations 
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Substituting  (6.2)  into  (6.1)  and  using  the  results 

_3_  _ k _3_  _3_  _3_  _3_ 

3r  m 3yx  ’ 30  = t 3s  “ m 3£ 


-^G(s)u(£) 

r=R  k s 
t 


we  obtain  the  magnetic  field  on  the  surface,  namely, 

Hr(Rt,e,4>)  ^ 0 


(6.2a) 


(6.2b) 


(6.2c) 


(6.3a) 


(6.3b) 


H (R  ,0,4>)  'v*  -sin 
v t 


♦I]1  - kt)v<°  - + “JG<S> 


where 


v«) 


(6.4c) 


(6.5) 


From  the  derivation  of  U and  V in  (6.1)  it  may  be  shown  that  (6.4)  is  asymp- 

-2/3 

totically  accurate  only  up  to  and  including  terms  of  order  (kR£)  . Thus, 

the  last  terms  in  (6.4b)  and  (6.4c),  which  are  of  order  (kRt)  may  or  may 
not  be  completed.  We  hereby  modify  the  factor  M in  (6.5)  to  read 


M - JS-  > 
kRt 


v’(0  + u’U) 


(6.6) 


The  only  "justification"  for  replacing  (6.5)  by  (6.6)  is  that  this  replace- 
ment yields  no  appreciable  difference  in  the  sphere  problem  (where  R /R  = 1), 

t b 

but  significantly  improves  the  numerical  solution  for  the  cylinder  problem 
(where  R.  /R^  may  be  very  large)  as  shown  in  Section  5.  We  have  not  yet 
succeeded  in  finding  a rigorous  justification  for  using  (6.6).  Referring  to 
the  cylinder  geometry  sketched  in  Fig.  1,  the  following  substitutions  are  made 

Hr  -*■  Hn,  (Hg/cos  <p)  -*•  Ht>  (H^/sin  $)  -*•  -Hb,  <f>  -*  6 . (6.7) 

Then  (6.4)  and  (6.6)  become  identical  to  (2.6). 

It  should  be  mentioned  that  expressions  similar  to  (6.4)  were  given  by 
Wait  in  1956  [15],  and  by  Hasserjian  and  Ishimaru  in  1962  [14],  [16].  Their 
expressions  contain  only  (ks)  ^ terms,  and  therefore  only  the  "hard"  function 
v(0  is  used. 


For  a given  magnetic  dipole  located  on  a perfectly  conducting  cylinder, 
the  surface  magnetic  field  is  given  in  (2.6),  which  is  an  approximate  asymp- 
totic solution  valid  for  large  kR,  and  may  be  used  for  any  observation  point 
on  the  cylinder.  The  derivation  of  (2.6)  is  based  on  the  following 
observation:  the  dominant  contribution  of  the  surface  field  depends  on  the 

curvature  of  the  surface  in  the  longitudinal  direction  of  the  ray,  not  that 
in  the  transverse  direction.  Hence,  we  simply  adopt  the  classical  solution 
of  Fock  for  a dipole  on  a sphere  for  the  present  cylinder  problem.  To 
include  the  effect  of  the  curvature  in  the  binormal  direction  of  the  surface 
ray,  we  have  arbitrarily  replaced  a factor  M of  the  Fock  solution  in  (6.5) 
by  that  in  (6.6).  A rigorous  justification  of  this  replacement  is  yet  to 
be  found.  A remarkable  feature  of  the  solution  in  (2.6)  is  that  in  the 
limit  kR  -+  «,  it  becomes  identical  to  the  known  exact  solution  of  a dipole 
on  a flat  ground  plane.  The  application  of  (2.6)  to  the  mutual  admittance 
calculation  yields  excellent  numerical  results,  and  therefore  (2.6)  may  be 
regarded  as  an  improvement  over  two  previous  asymptotic  solutions.  A 
future  research  problem  is  to  generalize  (2.6),  according  to  the  recipe  of 
GTD,  so  that  it  may  be  used  for  an  arbitrary  convex  conducting  surface. 
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APPENDIX 
FOCK  FUNCTIONS 

In  this  appendix  we  define  and  list  some  useful  formulas  of  the  functions 
w^(t),  w2(t),  v(£)»  u(£),  and  v^(C).  These  functions  are  commonly  known  as 
Fock  functions. 

(i)  Definition:  For  a complex  t and  a real  £, 


w.(t)  - — [ dz  exp  (tz  - 7:  z3) 

1 /iT  J r 


(A-l) 


w2(t)  = - 
/tt 


(1  3 \ * 

tz  - — z 1 = W^(t) 


jr/41/2  1 

’ w2(t) 

- e t,  — 

/n  ■ 

r w'(t) 
1 

j3n/4  3/2  1 

f W2(t) 

; 5 r,  J 

L w,(t> 
1 

j 3ir/4  3/2  _1_ 

w2( 

^ J 

rL  ' "F 

and  r2  are  sketched 

e'^  dt 


(A- 2) 


(A-3) 


(A-4) 


(A-5) 


derivative  of  w2(t). 

(ii)  Zeros  of  w2(t)  and  w2(t):  They  are  given  by 

1 - S,  * Itni6"^73  ’ and  t - t;  - , 

respectively.  The  magnitudes  of  the  first  ten  zeros  are 


(A-6) 


— 
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» lt„l 

i*;i 

1 2.33811 

1.01879 

2 4.08795 

3.24820 

3 5.52056 

4.82010 

4 6.78671 

6.16331 

5 7.99417 

7.37218 

n 

!‘„l 

i*;i 

6 

9.02265 

8.48849 

7 

10.04017 

9.53545 

8 

11.00852 

10.52766 

9 

11.93602 

11.47506 

10 

12.82878 

12.38479 

(iii)  Residue  series  representation:  For  real  positive  £, 


rjn/4A-  £1/2  i (t’)_1e 


(A- 7) 


vx(0 


2jW42^  c3/2  £ 

n=l 

>/42^  £3/2  £ 

n~l 


(A- 8) 


(A-9) 


v'(0  - | e"jl,/4/ir  r1/2  I (1  - j2St;)(t;)-1e  n 

n=l 


(A- 10) 


u'<0  3 ej7T/43^  £1/2  I (l  - j | 


(A-ll) 


(iv)  Small  argument  asymptotic  expansion:  For  real  positive  £.  and  £ -*■  0, 


»({)  4.  1 - & 'i'Hjn  + 11  {3  + e-J./4t9/2  . 4 U1  , 10-3{6  + _ 

bU  (A-12) 

te\  n i ^ jw/4  3/2  . 51  .3  . 5/n  -jir/4_9/2  _ n-2  6 

u(£)  ^l--T-eJ  £ + r±-  £ + T7~  e J ' £ ' - 3. 701  x 10  £ + . . . 

2 12  64  (A-13) 


,1<E)  4.  1 + f e^V'2  - ii  C3  - *5  e^'V'2  + 4.555 


x 10_2£6  + . . . 


v'<£)  * 
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2.485  x 10-2C5  + . . 


(A-14) 

(A-15) 
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u'(£)  - | ^ V/Z  e'^'V 


.221  * 10-1£5  + ... 

(A-16) 


(v)  Relation  to  Fock's  attenuation  factor:  The  attenuation  factor  V(x,y,q) 

is  defined  by  Fock  in  p.  207  of  [12],  and  it  reads  (after  replacing  i by  -j) 


, 1tt/4  1/2  1 f W2(t  " y)  -ixt  , 

‘’y’q)=e  X ^Jr  w*(T)  ■ -qw2-(t)  e dt 


(A-17) 


The  functions  v and  u are  related  to  V by 


v (x)  = 2 v (x,y  = o,  q = 0) 


(A-18) 


u(x)  = jx  lim  q — 

q-*30  L Z)  y=0 


(A-19) 


(vi)  Relation  to  functions  defined  by  Logan  [17]: 


v(£)  =2  5 nO 


(A-20) 


. . j3n/4 ,3/2-.  . 

u(£)  = eJ  £ n£)  • 


(A-21) 


(vii)  Tabulation:  Functions  w^t),  w2(t)  and  their  derivatives  are  tabu- 

lated in  [12],  while  i|j(£)  and  ij i(£)  are  tabulated  in  [17],  Numerical  curves 
of  u,  v,  and  v^  can  be  found  in  [4],  [14],  [15]. 

(viii)  Numerical  evaluation:  For  £ _>  £q,  the  residue  series  representation 

with  the  first  ten  terms  in  the  summation  may  be  used.  For  £ ^ £q,  the  small 
argument  asymptotic  expansion  with  the  first  five  terms  may  be  used.  It  has 


been  indicated  in  [4]  that  the  smoothest  crossover  is  obtained  if  £^  = 0.6. 
In  the  present  study,  we  set  £^  = 0.7,  where  the  difference  in  the  two 
representations  is  as  follows. 
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